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Abstract 

Motivated by the fully gapped superconductivity in iron-based superconductors with uncom- 
pensated electron pockets, we propose a spin-singlet analogue of the superfluid phase of 3 He-B. 
We show that the staggered tetrahedral structure of these systems with 2 iron-atoms per unit cell 
introduces additional site and orbital band indices r and a which act as internal isospin degrees 
of freedom, opening up a new class of tetrahedral and orbital ("TAO") triplet condensates. This 
permits two locally d-wave pairings to combine in quadrature without breaking lattice or time 
reversal symmetry, to produce a fully gapped topological superconductor. The internal chirality 
of this state is selected by a Josephson coupling between the d-wave components. This state is 
predicted to give rise to Raman-active excitations associated with the reversal of the gap chirality. 
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The discovery of superconductivity with T c = 26K in LaFeAsO by Hosono et. a/.[l| has 
opened a new field of iron (Fe)-based multi-band high temperature superconductors (SCs)j^, 
31. Transport and ARPES experiments consistently show that these systems are fully 
gapped, while optica, m ea_ts a ls0 snow a signet Coulo mb i nt e r ac ti o n fl. Mo, 
electronically mediated superconductors, avoid the onsite Coulomb interactions by building 
nodes into the pair wave, so that we expect the onsite pairing (^^)p e = J k <^ = 0. In the 
organic, heavy fermion, cuprate and strontium ruthenate superconductors, this "Coulomb 



orthogonalization" is guaranteed by a d- or p-wave symmetry of the order parameter 



Superconductivity in the Fe-based systems is widely believed to derive from a gap function 



with s^ 1 symmetry 2|,[l0|, where the gap takes the form A(k) = Aq + 2Ai cos k x cos k y , giving 



rise to a gap node between the electron and hole pockets. This is supported by experiments 



showing c-axis Josephson tunneling H|] and QPI measurements on Fe(Se,Te) 12|, including 



absence of the Wohlleben effect in SQUID measurements 13f] . In this scenario, a phase 



cancellation between the electron and hole pockets eliminates on-site pairing, minimizing 
;he Coulomb interaction. However, when the hole pockets are absent, such as in A x Fe2Se2 
and single layer FeSe 15], the node in the s ± order parameter is expected to intersect 
the the electron pockets to achieve Coulomb orthogonalization. These nodes have not been 
observed. Another possibility is that the symmetry of the order parameter changes in 
the electron pocket materials, for example by the development of d-wave pairing, with 



nodes lying between the pockets la]; yet this scenario appears to be inconsistent with the 



observation of small fully gapped hole pockets around the Z-point in Tlo.63K .37Fei. 78 Se2 



li 



These observations motivate us to seek an alternative description of the of the pairing 
symmetry which can account for the Coulomb orthogonalization in both the electron-hole 
and single electron-pocket materials without a change in condensate symmetry. In this 
search, we are inspired by superfluid 3 He-B, where the Fermi surface is fully gapped, yet 
the pair wavefunction eliminates local pairing through the formation of p-wave spin triplet 
pair condensate. In two-dimensional 3 He-B, the gap function takes the form (ck CT c_kcr'} = 
(k x (j x 4~ ky a y )ia 2 ; although the odd-parity k x and k y triplet components contain nodes, their 
anti-commuting spin structure causes them to add in quadrature to create a fully gapped 
condensate where A(k) oc [k% + ky]. While the spin-triplet scenario is ruled out in iron- 
based superconductors by Knight-shift measurements [la Q , the multi-band nature of the 



iron-based superconductors may offer an analogous route to a fully gapped condensate, for 
as we now show there are two additional quantum numbers, r = ±1 and a — ±1, arising 
from orbital and lattice symmetries, which open up a new class of non-trivial triplet states. 

The staggered tetrahedral structure of the iron based superconductors results in an en- 
larged unit cell containing two rotated Fe tetrahedra. The site index associated with this 
novel structure plays the role of an additional internal "tetrahedral" index r = ±1 labeling 
electrons on iron sites within "right" or "left" pointing tetrahedra (Fig.[TJ) The Pauli isospin 
linking these states will be denoted by f = {ti,t 2 ,t 3 ). Likewise, we associate a quantum 
number a = ±1 with the two degenerate zx (a = +1) and zy (a = — 1) d-orbitals of the 
iron sites, denoting the Pauli matrices that link these states by a = («i, a 2 , a 3 ) (Fig. [TJ). 




T = +l T = —\ a = +l a = -l 



FIG. 1. (a) Tetrahedral index r = ±1 denotes Fe in right or left pointing tetrahedra, (b) orbital 
index a = ±1 denotes zx and zy orbit als. 

The electron operators now carry three discrete quantum numbers, so we replace the conven- 
tional electron operator Cko- — > c^, where v = (r, a, a) denotes the triplet of site, orbital and 
spin indices. Since we are treating orbital triplets, we adopt a generalized Balian-Werthamer 
spinor to describe the electron fields 

where A2 = T20t2&2 flips the tetrahedral, orbital and spin quantum numbers so that the 
electron and hole fields transform the same way under the various isospin rotations. In this 
notation an extended s ± pair, written A s± (k) = cos{k x )cos{ky)T2Q.2^-ai is an orbital and 
site triplet, the orbital analog of equal spin triplet pairing, for once multiplied by A2, it is 
diagonal in site and orbital indices. We now consider a more general class of orbital triplet 
pairing. 
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The 2D Fe As(Se) layer exhibits C& v symmetry about the As(Se)atoms, but the two iron 
atoms in the unit cell lie at sites of lower symmetry. There is no 90° rotation symmetry 
about an iron site - instead, the unit cell is invariant under a combined 90° rotation and a 
translation between the two iron sites. This means that the site and orbital indices must 
transform non-trivially under symmetry operations of the C± v group. Moreover, the global 
symmetry of the superconducting state need not correspond to the nominal spatial symmetry 
of the pairs at the iron sites. Thus for example, one can have a staggered d-wave paired state 



which has a global s-wave symmetry [201] . Lastly, we note that there is no local inversion 
symmetry at each iron site; instead the inversion is located mid-way between neighboring 
Fe atoms. 

We shall argue that the strong Coulomb forces at the iron sites favor states with local 
(i-wave symmetry, seeking gap functions which contain terms with both local d xy and d x 2_ y 2 
spatial symmetries coupled to different site and orbital channels. Inspired by 3 He-B, the 
pair wave-function then takes the following form 

(c ki /C_ kl /) = (A x . y (k)r^„ + A x a_ y a(k)I^ w //)(iA 2 )^v, (2) 

where A xy (k) = Aj sin k x sin k y and A a ,2_ J/ 2(k) = A 2 (cos k x — cos k y ) are the xy and x 2 — y 2 
pair amplitudes and the matrices T 1 ' 2 £ [r a <8> cc&j are tensor products of the site and orbital 
Pauli matrices. These spin singlet pairs have odd parity under spin exchange, but are 
even under momentum inversion, site (r) and orbital (a) exchange, thus preserving the 
antisymmetry of the wavefunction. 

To preserve lattice symmetry, Ti and T 2 must belong to the same irreducible representa- 
tion of the D^h group and in addition, both must have the same parity under time-reversal 
and inversion symmetry. The most important requirement, is that the two d-wave pairing 
components of the Hamiltonian anti-commute with each other, which guarantees that the 
two gaps add in quadrature. We call this "tao pairing" (tetrahedral and orbital: raa). 
One of its features is that as in 3 He-B, the d-vector d(k) = (A^k), A !B 2_ Jf a(k))/|A(k)| has 
topological character, in this case, winding twice around the origin as k runs around the T 
point. The chirality v = — sgn(A!A 2 ) of the order parameter in internal tao space is selected 
by the underlying hamiltonian without breaking any symmetry. 
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Irred. Rep. 


SC State 


Local Symmetry 


Orbital (a) Symmetry 


Ai 


c x c y T 2 a 2 


8 ± 


11 + 22 


Bi 


iA x 2_ y 2 r 3 a 2 


d x 2 —y 2 


11 + 22 


B 2 


A xy r 2 a 2 


d X y 


11 + 22 


Ai 


A a ,2_ J/ 2r 3 ai 


d x 2 —y 2 


11 - 22 


iA xy r 2 a 3 


d X y 


12 + 21 


A 2 


A x 2_ y 2T 3 a 3 


d x 2 —y 2 


12 + 21 


iA xy r 2 ai 


d X y 


11 - 22 


A\ 


iA xyz T\a 2 


3D staggered d xy ("7?" ) 


11 + 22 


B{ 


A xyz na>i 


3D staggered d xy 


11 - 22 



TABLE I. Symmetry classification of the singlet superconducting order parameters in the Fe- 
based systems. The non-trivial A\ and A 2 representations contain two one-dimensional d-wave 
components that transform identically under the C^ v group. The A* and B* representations 
involve staggered d xy order. The t\ operator which implements the staggered pairing is odd- 
parity under inversions through the center of inversion, so that the "f-wave" combination A xyz r\ = 
(A xy Ti) x sin£; 2 is even parity. Since a 2 and t 2 are odd under time-reversal, the time-reversal 
invariant combinations ia 2 and ir 2 are used in the pair operator. 

Table. |T] classifies the even-parity order parameters, with nearest or next-nearest neighbor 
pairing in the Fe-based SCs according to the irreducible representations of the C± v group. 
Within the A\ and A 2 symmetry classes there are two one dimensional representations 
with local d xy and local d x 2_ y 2 symmetry. Their identical symmetries under C^ v means that 
hopping terms in the Hamiltonian will interconvert the two types of pair, creating an internal 
Josephson coupling between them. Remarkably, as we shall now show, the resulting presence 
of two linear combinations of <i-wave states with different local spatial symmetries, but the 
same global symmetries, i.e. the same irreducible representations, leads to the formation 
of a topologically non-trivial state in which the gaps add in quadrature, eliminating the 
underlying nodes from the quasiparticle excitation spectrum. 

In this paper, we shall focus on the A\ and A 2 states, which are invariant under inversion 
symmetry of a single layer. By contrast, the A \ and B{ representations, which involve 
staggered d xy pairing, are odd-parity within a single plane. These states can only form a 
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tao state in bulk 3D systems, where the inversion center lies between FeAs(Se) layers; by 
staggering the d xy pairing between layers, the A\ and B\ states acquire an additional sin k z 
dependence, restoring the inversion symmetry and giving rise to the A xyz = (staggered A.,,,,) x 
sink z form-factor. We note that the A\ state corresponds to the //-pairing proposed in 20 ]. 
We shall return to disucss this case later in the paper. 

Using this scheme, we are able to identify two fully gapped 2D tao states, one with A\ 
symmetry, the other with A2 symmetry (see Table. |TT])- In multi-layered iron-based systems, 
the inversion symmetry of an isolated layer is broken, allowing us to also consider tao 
states which mix "starred" and "unstarred" representations with opposite in-plane inversion 
symmetry. 



tao state 


Irreducible Representation 


Description 


A^r 2 a372 + A x 2_ y 2T 3 aiji 


Ai 


Uniform orbital triplet. 




A 2 


Uniform orbital triplet. 


^xyzTia 2 -j2 + A rr 2„j / 2T3ai7i 


A\ 


Staggered d xyz orbital triplet. 


A-xyznaill + A x 2„j / 2T 3 a272 


B{ 


Staggered d xyz orbital triplet. 



TABLE II. Family of orbital triplet tao superconducting states. The asterix denotes tao-states 
which break the inversion symmetry of a single FeAs layer. The Nambu matrices acting in particle- 
hole space are denoted by 7 = (71, 72, 73). 

To illustrate these results, we now carry out a BCS treatment of the Ai tao paired state 
(see Table. [II]), with BCS model 



(3) 



k k,k' 

where N q is the number of sites in the lattice and 



k,k' 



bl 



kl 



c 



d xy (k)iT 2 a 3 )i\ 



c 



-ki/'' 



bl 



k2 



"ki/ 



4.2^2 (k)r 3 ai Iz'A 



-ki/' 



(4) 



create the two components of a tao-pair with A\ symmetry, where d xy (k) and d x 2_ y 2 (k) are 
the G?-wave form factors. For simplicity, we have omitted interaction terms that mix the the 
two pair types, since as we shall see, this mixing is already provided by the hopping. 
The kinetic part of the Hamiltonian 



h(k) = e (k)7 3 + ei(k)n + e 3 (k)«i + e 4 (k)r 1 a 3 73, 
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(5) 



contains the kinetic hopping terms consistent with the symmetries of a single FeAs layer. 
Denoting q = cos A;/ (/ = x, y) and si = sin ki (I = x,y), then ei(k) = 4ti(c x + c y ) and 
eo(k) = AtQC x Cy — fi are the intra-orbital (orbitally identical) nearest and next nearest- 
neighbor hopping, where \i is the chemical potential. The terms e^k) = At±{c x — c y ) and 
e 3 (k) = 4tsS x s y are the inter-orbital nearest and next-nearest neighbor hoppings (we shall 
choose a gauge in which t^t^ > 0). These terms act as tetrahedral and orbital fields that split 
the orbital degeneracy of the electron and hole pockets. The resulting normal state spectrum 
is given by e nQ (k) = e n (k)+a v / e 3 (k) 2 + e 4 (k) 2 , where e n (k) = e (k)+sgn(n)ei(k) and n, a = 
±1 are two sets of band indices. Since e 3 (k) = &t$d xy ()s) and e 4 (k) = At^d^ _ y 2(k) have d- 
wave symmetry, the product of the these terms generates matrix elements that hybridize 
the two locally d-wave components of an A\ or A2 tao pair. 

Performing a mean-field decoupling of the interaction, we obtain 



V \9192J 



(6) 



where 

H(k) = h(k) + A :r j / (k)r 2 a 3 72 + A a ,2_j /2 (k)r 3 ai7i (7) 

is the Nambu Hamiltonian. "H(k) can be diagonalized, giving rise to four separate quasi- 
particle eigenvalues concentrated around the original electron and hole pockets, given by 



rpna 



A n (k) + a v / ^n(k) 2 - 7 n(k) 2 



1/2 



(n,a = ±l) (8) 



where 



An = e 2 n + el + e 2 4 + A 2 xy + A 2 x2 _ y2 , (9) 
ll = (4 -4-4 + + A' 2 -, 2 ) 2 + 4(A xy e 3 - A x2 „ y2 e 4 ) 2 , 

and for clarity we have suppressed the explicit momentum labels k. From (jSJ) and (Q we 
see that (E£ + E£~) 2 = 7„,(k) 2 > is positive definite, so that the excitation spectrum is 
fully gapped. In fact, there are two types of gap structure, depending on the chirality 
v = — sgn(AiA2) of the gap: for v < 0, (Ai, A2 in phase), the second term in 7 n vanishes 
on the Fermi surface and the gap becomes highly anisotropic, almost closing with a very 
small minimum value A min ~ min(Ai, A 2 ) 2 /W, where W is the bandwidth; for positive 
chirality v > (Ai, A 2 out-of-phase) , all terms in 7 2 remain positive, and the gap is weakly 



anisotropic, with a large minimum value A m ; n ~ min(Ai,A2). The large gap, positive 
chirality state is thus selected by the lattice. Had we not chosen a gauge where £3^4 > the 
selected chirality would be given by v = sgn(t 3 i 4 ). 

Figs. |2] & [3] show the superconducting gaps for the two generic cases seen in the iron- 
based superconductors; two electron and hole pockets around M and T, and also two electron 
pockets around M, showing the fully gapped structure in both cases. To explore the internal 



a) b ) 
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FIG. 2. (a) Superconducting gap of the hole and electron pockets around T and M. The Fermi 
surface contours are shown as red lines, and the superconducting gaps are shown in red and blue 
for the two hole pockets, and in yellow and green for the electron pockets. (b)The DOS for the 
superconducting state with a clear gap around Ep. Energies are measured in units of to- The 
parameters used are to = 1-0, ti = 0.1, £3 = 1.0, £4 = 0.6, /i = 0.6, Ai = 0.5, A2 = —0.2. 



Josephson coupling between the two condensates, we minimize the mean-field Free energy 



F = N„ 



A? A 2 2 
— + — 

9i 92 



2T J> 

k,n,a 



2 cosh 



2T 



with respect to Ai and A 2 , to obtain two coupled gap equations 

1/gi-Xxy -x.j 

-X.J 1/92-Xxt-y 2 



)l 


A1 l 




\A 2 



where, denoting J k = / ^ and B kn = A\ n - 7^, 



th(/3E£72) 



1 + 2a 



(4t 3 d x 2 



B 



(10) 



(11) 



(12) 



S 



FIG. 3. (a) Superconducting gap for a choice of parameters with only electron pockets around M. 
The Fermi surface contours are shown as red lines, and the superconducting gaps are shown in 
yellow and green for the two electron pockets, and the black lines mark the real Brillouin zone. 
(b)The DOS for the superconducting state with a large gap around Ep. Energies are measured 
in units of to- The parameters used are to = 1.0, t± = 0.15, t$ = 0.7, t$ = 0.7, fi = —4.0, Ax = 
0.5, A 2 = -0.25. 



Xx 2 - 



th(/3E^/2) 

2E na 



2a 



(4t 4 d. 



j-yj 



B 



kn 



J 2 



are the d xy and d x 2_ y 2 pair-susceptibilities, while 



Xj 



E 



a 



th(/3E°"/2) 



2Ut 



3 t 4 



B 



kn 



[Adxyd^ 



(13) 



(14) 



is the internal Josephson coupling between them. Since the E^ a have the lowest magnitude 
when a = —1, xj is negative, inducing a ir coupling between the d-wave condensates, 
favoring selection of a positive chirality v > 0. The mean-field transition temperature is 
determined by (1/gi— Xxy)(]-1 9i— Xx 2 - y 2 )— Xj = 0. In a conventional d-wave superconductor, 
pairing develops exclusively in the single most attractive d-wave channel, but in this system, 
the strong Josephson coupling between the two states will immediately lock the two gaps 
into a topological state with a single, enhanced transition temperature. To get an idea of 
the magnitude of this effect, let us put Xx y ,Xx 2 - y 2 ~ N(0)ln(W/T), Xj(T) ~ — \jx xy {T), 
where iV(0) is the density of states and Aj a quantity of order unity, then if g x = g% = g, 
the transition temperature is given by T c ~ W exp [— 7 ) g N(o) 1 ' so ^ ne m t erna l Josephson 
coupling enhances the coupling constant gN(0) by a factor (1 + Aj). 



In a strong coupling picture, tao-pairing is expected to be driven by superexchange pro- 
cesses, giving rise to Kugel-Khomskii type J\ — Ji interactions involving both spin and orbital 
exchange 21]. A decomposition of these terms will generally give BCS pairing terms in all 
the irreducible representations of the point group symmetry of the lattice, but the A\ and 
A 2 tao states are then expected to exhibit the highest temperature instability thanks to the 
enhancement effect of the internal Josephson coupling. 

We end with a discussion of the possible implications of tao-pairing. The tao scenario 
provides a mechanism for understanding how singlet superconductivity can avoid the strong 
Coulomb repulsion on the iron sites, namely through the locally (i-wave character of its two 
pairing channels. At the same time, the scenario naturally accounts for absence of nodes in 
the 11-systems that only have electron pockets around the M points. There are a number 



of well established features o: 
peak in NMR measurements 



the experiments, such as the the absence of a Hebel-Slichter 



scattering 



22 



19| , and the observation of resonance peaks in the neutron 



23] that are qualitatively consistent with the ci-wave character of tao pairing, 



but have yet to be studied in detailed calculation. We have also not yet carried out a detailed 
study of the interaction of tao pairing with the Coulomb interaction. 

In 3D multi-layer iron systems, tao pairing can in principle occur in four different global 
symmetry representations, A±, A 2 and the staggered A\ and B\ representations, driven by 
interlayer physics that break the local inversion symmetry of the layers. The uniform A g 
(s-wave) symmetries would appear to be the best candidate scenarios: driven by a strong 
in plane Josephson coupling, rather than a weaker inter-plane effect we expect these paired 
states to exhibit the greatest enhancement of T c . Moreover, this scenario allows us to 
understand the observed c-axis tunneling into s-wave superconductors fll|. The A\ state 
appears to possess the most isotropic gap (see Fig. m&E]), and would appear most consistent 
with STM experiments 24]. 

The staggered d-wave A\ and B\ states may become relevant to systems with strong 
interplane coupling. Many ex per iments suggest a multi-gap character to the superconducting 
condensate, including STM 



24. 



and the observation of a 



19] 



25|, ARPES experiments 
cross-over in the NMR relaxation rate from from T 3 at high temperatures to T 5 
behavior at low temperatures. This is consistent with the formation of a staggered A\ or 
B{ tao condensate in which the presence of the staggered d-wave gap A xyz = Ais x s y s z will 
lead to out-of-plane point nodes. The development of these nodes would naturally account 
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for the observation of a T 5 dependence in ^ measurements in NMR experiments. 

One of the interesting consequences of the theory, is that it predicts the condensation 
of pairs with an internal chirality v. Pairs of reversed chirality, generated by anti-phase 



fluctuations of the two gaps are expected to then lead to a low-lying Leggett mode [27|, [28] 
that will be Raman-active. There should also be higher-energy excitations of Cooper pairs 
of a different symmetry, for e.g. from A\ to B\, which will couple to corresponding Raman 
modes. The emergence of these low-lying "failed d-wave" excitations below T c is key property 
of the tao-condensate. 

A second aspect of the hidden chirality, is its topological character. An analysis of the 
symmetries of the pairing Hamiltonian shows that it lies in the "CI" class of topological 



superconductors 



29] . In two dimensions, these states are topologically trivial and are not 



expected to exhibit edge states. However, the three dimensional extension of these pairing 
states does contain a topological Z invariant which allows for the possibility of non-trivial 
surface Andreev states. 

In conclusion, symmetry analyses of the iron-based superconductors lead us to predict 
the existence of a new class of spin singlet, tetrahedral and orbital triplet pairing ("tao 
pairing"), in which a topologically non-trivial combination of locally d xy and d x 2_ y 2 paired 
states gives rise to a fully gapped excitation spectrum. This novel state utilizes two new 
quantum iso-spins, providing an interesting new candidate symmetry for the fully gapped 
iron-based superconductors. 
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